Abstract 



We derive a determinant formula for the WKB exponential of singularly perturbed 
Zakharov-Shabat system that corresponds to the semiclassical (zero dispersion) limit of 
the focusing Nonlinear Schrodinger equation. The derivation is based on the Riemann- 
Hilbert Problem (RHP) representation of the WKB exponential. We also prove its 
independence of the branchpoints of the corresponding hyperelliptic surface assuming 
that the modulation equations are satisfied. 
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1 Introduction 

The semiclassical analysis of the focusing Nonlinear Schrodinger (NLS) 

ied t q + ^e 2 d 2 x q+\q\ 2 q = 0, (1) 

has produced [21 [3j H] solutions that one recognizes as modulated multi-phase periodic or 
quasi-periodic waves. These wave solutions of NLS are expressed in terms of hyperelliptic 
theta functions (see PQ) built from the radical 

1/2 



R(z) 



'4JV+1 

Y[(z- on) 

. i=0 



(2) 



where even a 2 k he in the upper halfplane and a 2 k+i = «2fc- Let 1Z = 1Z(x,t) denote the 
hyperelliptic Riemann surface of R, where oriented arcs 7^ fc , connecting and a^, 

k = 1, 2, • • • ,N, their complex conjugates 7~ k , connecting a^+i and a^k-i, together with 
jmfl connecting a± and a>o, form the branchcuts (main arcs), see Fig. [U 

Points on depend on x, t but do not depend on e. They are called branchpoints of the 
hyperelliptic Riemann surface TZ or simply branchpoints of R. The number of wave-phases of 
a solution q(x, t, e) of ([I]) (in the limit e — * 0) is equal to the genus 2N of the corresponding 
Riemann surface lZ(x,t). The branchpoints satisfy a system of equations known as the 
modulation equations or the modulation system that is discussed below. 

Nonlinear Schrodinger equation is one of the most celebrated examples of an integrable 
PDE, i.e., a nonlinear PDE that can be "linearized" through the Lax pair. The u x" (spatial) 
part of the Lax pair is a second order system of linear ODEs 

*,* = --(-, \ , ^ £ ))<S>. ( 3 ) 
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q(x,t,e) —z 
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Figure 1: Riemann surface TZ(x,t) 



known as Zakharov-Shabat system (see [5]). In the limit s — > (semiclassical limit of the 
NLS) system ([3]) becomes a singularly perturbed system, and, as such, is a subject of the 
WKB analysis. 

The goal of this paper is to find a determinant formula (see ffTB|) ) for the WKB exponen- 
tial g of to prove its independence of the branchpoints a assuming that a satisfy the 
modulation equations, and to derive various forms of modulation equations for a using the 
determinant formula. Our derivation is based on the RHP representation (jSJ) of g, which 
was obtained and discussed in [I] for (jTJ) with pure radiational or radiational and solitons 
initial data (for the pure soliton case, see [3J). Connection between the WKB exponential g 
and the RHP (j5j) is studied in a separate paper [?]. The RHP (j^J) and its solution are the 
main objects of this paper. 

The inverse scattering method of integration of (jTJ) is based on the scattering transform, 
i.e., on the connection between the initial data of ([1]) and the scattering data of ([3]). In 
general, the scattering data consists of the reflection coefficient defined on K. and of the 
eigenvalues of Zakharov-Shabat system together with their norming constants (these eigen- 
values correspond to solitons). In the semiclassical limit of (JTJ), the reflection coefficient 
r (z,e) depends on z and e. We denote 

r(z,e) = r (z,e)e^ (xz+2tz2) and f = -eihxr. (4) 

In general, / is a function of z and e. However, studying the semiclassical limit of (jTJ), we can 
consider only the leading order term (in e) of /, see [I]. Throughout this paper, we assume 
that f(z) has analytic continuation into the upper halfplane with the exception of a finite 
number of logarithmic branchcuts and of isolated singularities. We further assume that the 
main arcs (branchcuts) of the hyperelliptic surface TZ(x, t) do not intersect the singularities 
of f(z). The values of f(z) in the lower halfplane are obtained by Schwarz reflection. Thus, 
in general, ^sf(z) has a jump on the real axis. 

To define modulation equations, we first assume that the branchpoints «2j — &2j(%,t), 
j = 0, 1, • • • , 2N, are known. Let r y ck be oriented arcs connecting «4fc-2 and a^k-i, a4fc-3 

respectively. These arcs are called complementary arcs. Let •y m ^ = J^k U 7m.fc? lc,k = 
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Figure 2: Contours 7 m ,2, %, 1 



1c k U 7 C fc> ^ = 1) 2, • • • , iV, where orientation is inherited. Define g(z) as the solution of the 
RHP: 



9+ + 9- = f + Wj on the main arc 7 m j, j = 0, 1, • - - , N 
9+ — g~ = Qj on the complementary arc j c j, j — 1, • • • , N 
g(z) is analytic at z = oo , 



(5) 



where all Wj, Qj are some real constants. We assume that the real constant Wo = on the 
main arc j m fi that connects ol\ and olq. 

Let 7 denote the unioun of all the main arcs 7m, j, j = 0, 1, • • • ,N and all the comple- 
mentary arcs with the inhereted orientation. It is well known that solution to the RHP §5§ 
is given by 



9{z) 



2ni 



no 



(C - z)R(0- 



1=1 \ J lm,i 



(C-z)R(C) + C + J JC:i ((-z)R(() 



dc, 



(6) 

Expressing the integrals over the arcs as integrals over the loops shown in Fig. [2j we obtain 



9W 



R{z) 

Ani 



Wi 



(C - z)R(C) 



,< (C-z)R(C) 



d( 



, (7) 



where the loops 7 and 7^ and the contours %,i (7^ consists of the sum of two arcs oriented 
oppositely as in the figure) are contractible to their corresponding arcs without passing 
through z. 

Deforming 7 so that now z is inside the loop 7 and still outside the loops j mji and 7 Cji , 
we obtain 



h(z) 



R(z) 

2m 



w t 

„, (C - z)R(C) 



7c, i 



(C - z)R(C) 



d( 



• (8) 
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where 

h(z) = 2g(z)-f(z). (9) 

The function h(z) is obtained by multiplying g by a factor of 2 and the residue — / being 
picked up as z cuts through the loop 7. 

If z approaches the ith main arc from either side, a residue is generated as z cuts through 
the loop 7 TO> j encircling the arc; multiplied by the factor outside the integral, the residue 
yields the contribution Wi to h. Similarly, if z approaches the ith complementary arc, the 
contribution to h from z cutting the contour 7 ci is +f2j or — f^, depending on whether z 
approaches from the left or right of the contour. These observations lead directly to the 
jump conditions 

h + + h- — 2Wj on the main arc j m j, j — 0, 1, • • • ,N 

h + — h- = 2Qj on the complementary arc j c j, j = 1, • • • , N 

h + — h_ = —2iQf on the real axis. (10) 

To see this, one takes into account the above residue calculations and the fact that the 
expression in the square brackets in (jHJ) is analytic in a neighborhood of any ctj (we are 
assuming distinct dj). From the above, it is clear that at any a = a 2 j 

h(z) ~ W + ±Sl + vi(z - a) 3 + u 3 (z - a)t H . (11) 

where ± refers to whether z is left or right of the contour and W, Q denote real constants on 
the main and the complementary arcs, adjacent to a = a^j. 

According to ([7j), g(z) ~ 0(z 2N ) as z — >• 00. The requirement that g(z) is analytic at 
z = 00 defines the system of 2N linear equations for Wj, flj 




fc = 0,l,--- ,2N-1. (12) 



Modulation equations comes from the requirement that the L 2 -solution to the RHP (jSJ), 
i.e., the Cauchy operator in (jH]), commutes with the differentiation in z. The equivalent 
statement is that at every a = a>2j the coefficient v\ = 0, so that ffTT]) becomes 

h(z) ~ Wj ± Qj + 1/3(2 - a)i H . as 2 -> a. (13) 

Equation f|T3l implies that the expression in the square brackets in equation (jSJ), let us call 
it B(z), vanishes at every «2j, i-6. 

(C - a 2j )/2(C) ^ \A m> < ^ ~ «2i)i?(C) 7 Tcji (C - a 2 j)^(C) / 

(14) 

where j = 0, 1, • • • , 2N and c^- is inside the loops around the main and the complementary 
arcs that are adjacent to a^j but outsisde all other loops and 7^. The modulation 
equations (fl4l is a system of 2N + 1 complex conditions satisfied by the 2N + 1 complex 
branchpoints a 2 j, j = 0, 1, • • • , 2N. It is also satisfied by their complex conjugates a 2 j+i- 
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2 Determinant formula 



To simplify the notations, we consider below the case N = 1. The obtained formulae allow 
a straightforward generalization to the arbitrary iVeN case. Let 



D 



r _d(_ r (d(_ 

J*f m R(Q Jim R(C) 

r jK_ r Cd£_ 

H R(0 ?% R(0 



K{z) = ±, 
2m 



7c R(0 

7 R(0 J 7 R(0 



j C/(CK 



dC 



$y™ (C-^R(C) 
^7. 

^7 (i.-.)/.'lU 



(C-*)fl(0 
/(CR 



(15) 



where 7 m = 7 m i and 7 C = 7 c x- It is well known that D ^ if all a 2 j, J = 0, 1, 2 are distinct. 



Multiplying the first two rows of K(z) by W = W — 1 and ft 
them to the third row and utilizing ([8]) and (|T2l) . we obtain 



f2i respectively, adding 



(16) 



where z is inside the loop 7 but outside the loops j m , %. That will be our standard assump- 
tion about the location of z for the rest of the paper, unless specified otherwise. It is clear 
that moving z inside the loops 7 m ,7 c would generate residue terms W and ±Q (depending 
on the direction z crosses the oriented loop %) in the right hand side of ffTB"]) . Combining 
this fact with ffT3l) . we obtain a new form of modulation equations 



K{a 2j ) = 0, j = 0,1,2. 
Lemma 2.1. Let a denote one of the branchpoints. Then 



dK{z) _ h(z) 



da 
1 

2wi 



D 



2(z 



a 1 



dD 

da 



;i7) 



(18) 



Proof. Let us write K(z) 
determinant K(z), see f|T5|) . Using 

d 1 1 



^1,^25-^3(^)1, where Kj denote the jth column of the 

c 



we obtain 
dKx 



da 



,K 2 ,K 3 (z) 



daR(C) 2((-a)R(() 



-\K 3 (a),K 2} K 3 (z) 



and 



C - a 



1 + 



a 



(-a ' 



(19) 



<9if 2 
9a 



- \Kx,aK 3 (a),K 3 (z) 



(20) 

Multiplying the first two rows of each of the above two determinants by W and Q respectively, 
adding them to the third row and utilizing ([HD, ffT2"l) and f|T7|) with a 2 j = we obtain (as 
before) 

Hz) 



1 

2wi 
1 

2?ri 



da 



K 2 ,K 3 (z) 



^,K 3 (z) 



(C-a)K(C) 



f 

$y e (C-a)fl(C) 
£ J% £ 

J^m R(C) J^r, 



R(C) 



7 

£y e RJO 



R(C) 



'7i 



(C - C T (C) 

7c (C-«)^(C) 



J2(*) 

J2(z) 



(21) 
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Adding these two determinants while taking into account (jT5l) yields 



1 

2tH 



8K X 



da 



K 2 ,K 3 (z) 



+ 



dD 
da 



Kz) 

R(z) 



Notice that 



follows from 
obtain 



dK 3 (z) 
da 



[K 3 (z) - K 3 (a)} 



(C-cO(C-z) 



1 

2~Tti 



K U K 2 , 



2(z-a) 

and (TTP1) . Using again (fT7j) with a 2 j 

dK 3 (z) 



1 

C-a 



2(z 



a 



-K(z) 



da 

Adding (I22p and (I24p completes the proof. 

Lemma 12.11 is the basis for the following theorem. 
Theorem 2.2. Let a denote one of the branchpoints. Then the equation K(a) 

d 



da 

Proof. According to ffTB"]) . we have 

dK(z) _ dh(z) 
da 

Substituting (HBJ) and (DSD into 



h(z) = 



D 



da R(z 
we obtain 
D d 



+ h(z) 



d D 

da R(z) 



R(z) da 



h(z) = 



which implies ( 1251) . 

Corollary 2.3. Modulation equations K(a2j) 

±h(z) = ^h(z), 



0,J 
d 



1, 2 imply 

d 



dt k[z) S di h{Z 



(22) 

(23) 
a, we 

(24) 
□ 

implies 
(25) 

(26) 

(27) 
□ 

(28) 



Remark 2.4. All the results of this section, unless mentioned otherwise, remain true for 
arbitrary genus 2N, N G N. Moreover, they do not depend on Schwarz symmetry of 7 and 
of f(z), as well as on any particular form of the functional dependence of f(z) = f(z; 0) on 
the external parameter(s) j3 (in the discussion above, (3 — x,t) and, in fact, are true in much 
more general setting. Indeed, let 7 = j(/3) be a Jordan curve in C and let f(z) = f(z; (3) be 
analytic (in z) on some open set S D 7 and smooth in j3. The contour 7 is partitioned into a 
finite number of interlaced nondegenerate (positive measure) main and complementary arcs 



by the branchpoints aj £ 7, j 

n 



1, 2, • • ■ , 2n. Then we have n main arcs j m j and n or 



1, depending on whether 7 is closed or not, complementary arcs 7 m ,j- 
the hyperelliptic Riemann surface 1Z((3) of the radical R(z) = \jYlj=i( z ~ a j 
g(z) = g(z] (3) satisfies the conditions (jSj), where Wj, flj are some complex constants. In fact 



The genus of 
is n — 1. Let 



all except any jV — 1 of these constants can be choosen arbitrarily. Let h(z) 
Then the modulation equation ( fT3l) at any a = aj implies = 0. 



2g(z)-f(z) 
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Since D and R do not explicitly depend on x, t, Corollary 12.31 together with (|T6l) imply 



that 

Using (TO)]) and 
where fo(z) 



±h(z) = ^—K(z) ±h(z) = 

dx { ) D dx [ h dt [Z) D dt 



K{z) 



/(*) = /„(*) 
~ei\nro(z), we calculate 



xz 



dt 



2tz 2 



£ * 

r dj 



(29) 
(30) 

(31) 



and 





dt 



K(z) 



y m ((-z)R(O 

; £ 

7c (C-z)R(() 



ML 

7c R(C) 



a 



j=0 



(C-*)«(0 

7c (C-z)R(O 



f 
f 



<-§EL °i<*C 



«(0 

K(C) 



(32) 

where z is outside the loops j m , %. 

It is easy to see that if z is inside loops 7 m ,7 c the j m integrals in (}3~Tj) . (1321 are equal 
to the corresponding integrals on the segment [o , o ] multiplied by —2, and % integrals in 
(I3"T]) . (1521 are equal to the corresponding integrals on the segment [0:4,0:4] multiplied by 2. 
If z is outside any of the loops j m ,%, then the corresponding residues should be taken into 
account. 



3 Differential form of modulation equations 

Modulation equations (j!7j) can be rewritten as ODEs 



dK 
da 



o 3 



d_ 

dx 



dK 
da 



a t 



dt 



(33) 



where a denotes the vector a = (ao, 02,04) (alternatively, we can consider a to be a 6- 
dimensional vector). According to (1T51) . in the case j 7^ I, j, I = 0, 2, 4, we have 



dKjon] 
da,- 



Hz) 



lim 

z^ai R{z) 



D 



2(z — a^ 



+ 



dD 

da n - 



(34) 



where z is inside any of the loops j m , 7 C that surround a m . Then \im z ^ ai jjA = 0. So, 
dK ( a ^ = 0. That means that the matrix is diagonal. 



To calculate dK( - a ^ we notice that 



1 

2tH 



da, ' 



daj 



, K 2 , K 3 (aj 



is dK * is ( - 



dD h(z) 
— — ■ hm 

daj z^oij R[z) 



. 



(35) 
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Let us take j = 2. Then 
dK(a 2 ) 



3D 



1 

/(C) 



-d( + 



dK 3 (a 2 ) 



W 



(C - a 2 fR{Q 

Equation f[T3"j) in a vicinity of z = a 2 can be rewritten as 

h(z) = W ± tt + c 2 (z - a 2 )R{z) + 0(z - a 2 



(C - a 2 ) 2 R(C) 



dc 



(36) 



(37) 



where c 2 G C. If z is inside any loops j m , %, the constants W, Q in ( 1371) should be replaced 
by zeroes, so we have 



R(z) 



c 2 (z - a 2 ) + 0(z - a 2 ) 2 . 



(38) 



Differentiating ( l38i) and taking into the account (JED, we obtain 

h(z) 



c 2 = lim 

2^«2 V-fl(2j 



I 

2txI 



f(0 



LJy 



(C - a 2 fR(Q 



W 



Thus, 



(C-a 2 ) 2 i?(C) 
dK{a 2 ) 



n 



7c 



(C - « 2 ) 2 i?(C) 



<9« 2 



-C 2 D . 

2 



(39) 



(40) 



Formulae (!36l) - P0l) are also applicable for a , a 4 , if formulae (|36|) . (I371) . (|39|) contain only 
integrals over the loops that are adjasent to a , a 4 respectively, and only the constants W, Q 
that correspond to these loops. 

According to (|36l) -( l40l) . (1331) can be written as 



3 d 
- Cj D(aj) x = -—K( aj ), 

Ordinary differential equations (14TI) imply 



Tt K{ai < 



{ait 



where, according to (|3~T|) . (1321) 
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f 
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(41) 



(42) 



f 
f 





(C- 


-a^)iJ(C) 


K- 


-aj)B(0 




(C- 




K- 


-a,-)fl(0 



Jy, R(C) 



J y m 
Jy 



(43) 



This is the Riemann invariant form of modualtion equations written as PDEs. 
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Alternatively, diffrential form of the modulation equations (|4"T1) can be obtained by dif- 
ferentiating (157)) and the corresponding equations at ao,a 2 - At a 2 we have 



A-h(z) = (W) x ± (fi), - ^ Cj R(z)(aj) x + 0(z - a 2 ), 

= (W) t ± (O), - -c^Xa^ + 0(z - a 2 ). 

Moving z inside any loops 7 m , j c that surround «2 will eliminate derivatives of W, Q in 
According to (l29j) . (151)) and (|32|) . these derivatives are 



(44) 



Wj; = — — 

D 



7c 



r(cy 



and 



47T? 



C 2 Sj=0 tt J 



7c 



R(0 



a 



2ttz 



47TZ 

IT 



7m 



C 2 Sj=0 a 3 

W) 



d(. 



Now equations (j4"Ij) follows from (jilj) and fl29|) . Equations (j4"5]) - (|4"6")l also imply 

8ir 2 



w x w t 



D 



(45) 



(46) 



(47) 



Finally, since the Cauchy operator for the RHP (J5]) commutes with differentiation, we 
have 



h\z) 



R{z) 
2m 



/'(C) 



(C - z)R(C) 



d(. 



(48) 



Combining this with (157)) yield h'(z) = [|cj + 0(y/z — aj] in a vicinity of z = ay. 
Thus, 



/'(C) 



37TZ L (( - Qtj)i2(C) 



Substitution of (HUD into flUD yields 



27r^^(a i ; 



no 



7 (C-ajWO 



in 



D§ f 



/'(C) 



7 (C-aj)fl(C) 



(49) 



(50) 



J =0,2,4. 
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